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ABSTRACT 

We present a model for the accretion flow around the supermassive black hole in the LINER nucleus 
of NGC 1097 which fits the optical to X-ray spectral energy distribution (SED). The X-ray segment 
of the SED is based on observations with the Chandra X-Ray Observatory , which are reported here 
for the first time. The inner part of the flow is modeled as a radiatively inefficient accretion flow 
(RIAF) and the outer part as a standard thin disk. The value of the transition radius (r tr ~ 225 Rs, 
where Rs = 2 GM/c 2 ) between the RIAF and outer thin disk was obtained from our previous fitting 
of the double-peaked Balrner emission line profile, which originates in the thin disk. The black hole 
mass was inferred from measurements of the stellar velocity dispersion in the host galaxy. When 
these parameters are used in the accretion flow model, the SED can be successfully reproduced, which 
shows that the line profile model and the accretion flow model are consistent with each other. A small 
remaining excess in the near-UV is accounted by the contribution of an obscured starburst located 
within 9 pc from the nucleus, as we reported in an earlier paper. The radio flux is consistent with 
synchrotron emission of a relativistic jet modeled by means of the internal shock scenario. In an 
appendix we also analyze the Chandra X-ray observations of the ~ 1 kpc circumnuclear star-forming 
ring and of an ultraluminous compact X-ray source located outside the ring. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: individual (NGC 1097) 
— galaxies: active — galaxies: nuclei 


1. INTRODUCTION 

Active galactic nuclei (hereafter AGNs) are thought 
to be powered by accretion onto a supermassive black 
hole. The accreting matter is thought to form a 
thin, equatorial accretion disk, which may reveal it¬ 
self observationally via its continuum emission (the 
“big blue bump”, produced by the sum of black body 
spectra fr om annuli in the disk at different tempera- 
tures; e.g.. lFrank. King fc Rainel2002HKoratkar fc Blaesl 
Another possible signature of a thin accretion disk 
is the profile of the Fe Ka X-ray line, which is emitted 
from the inner parts of the disk and distorted by strong 
relativistic effects ( see, for example, iFabian et al .11 19811 
iTanaka et, all 179951) . An alternative signature, broad, 
double-peaked Balrner emission lines from the outer ac¬ 
cretion disk (analogous to the lines of cataclysmic vari¬ 
ables) has been obse rved i n the optical spectra of some 
AGNs (see review bv lEracleouslll998lh 
Many AGNs with double-peaked Balrner lines (here- 
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after “double-peakers”) were disc overed in the spectro- 
scopic survey of radio-loud AGNs (lEracleous fc Haloernl 
M taxi ) an d by t he Sloan D i gital Sky Survey 
llStrateva et all 1200.11 iWang et all 1200511 . How¬ 

ever, NGC 1097 is an important object because 
it was the first (and the best studied yet) case 
of a low-luminosity LINER nucleus to display 
broad, dou ble-pe aked Ha an d H/3 emis sion lines 
llStorchi-Bergmann. Baldwin fc Wilsonl 11993 1. with a 
full width at half maximum of ss 7500 km s -1 . This 
discovery was followed by observations of such lines 
in other low-luminosity AGNs (LLAG Ns), mostly 


1996; iShields et all 120001: IHo et al.l 1200(1 iBarth et al 


with the Hubble Space Telescope ( HST, 

"all 


200 11. The appearance of the disklike emission lines of 


NGC 1097 is thought to be a transient event because 
they were not seen in previous observations. This 
transient event was possibly a tidal disruption of a 
star by a nuclear supermassive black hole, where the 
debris of the disrupted star feeds the accretion flow, 
although other possibilities were also discussed by 
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IStorchi-Bergmann et al.l (Il997f) . Since NGC1097 is a 
LLAGN, this transient phenomenon opened up the pos¬ 
sibility that we could witness variations in the disklike 
emission lines on a timescale of years or even months, 
and the variabi lity of the broad H a line s was indeed 
monitored llStorchi-Bergmaxm et a,1.1 ll 99,4 120031 here¬ 
after SB03). The observed Ha profiles of NGC 1097 were 
reproduced by kinemat ical m od els of non-axisymmetri c 
thin accretion disks bv IStorchi-Bergmann et all ( 1997 1. 
who inferred that the line-emitting disk is truncated at 
an inner radius of « 225 Rs, where Rs = 2 GM/c 2 is the 
Schwarzschild radius and M is the b lack hole mass, the 
black hole mass was measured by Iflewis Eracleousl 
(120051) as (1.2 ±0.2) x 10 8 Mg from the stellar velocity 
dispersion using the Ca II near-IR triplet lines. 

On the assumption that the broad double-peaked emis¬ 
sion lines observed in AGNs originate in a thin accre¬ 
tion disk, an external source of illumination is needed 
to power them. Some compelling arguments for this 
have been put forward by different authors, involv¬ 
ing the energy bud get of th e line- emitting region (e.g., 
iGhen. Halpern fc FilippenkdHl989i) and the low temper¬ 
atures achi eved b y the disk under local viscous dissipa¬ 
tion (lCollin-Souffrinlll987t) . For instance, considering the 
energy budget of the thin disk in a small sample of radio - 
loud double-peakers llEracleous fc I lalneriiltl 99l 12003 ). 
the disk can power its own lines only if virtually all the 
energy available locally is conver ted int o recom bi natio n 
l ines, which is unlikely (|Ghen. Halnern fe Filippenkoj 
H989|) . The presence of an inner, inflated, r adiatively 
ineffi c ient accretion flow (RIAF, see reviews byiQuataern 
120011 iNaravan. Mahadevan & Quataertlll998[ hereafter 
NMQ98) emitting an X-ray continuum which illuminates 
the outer thin disk solves all these energetic probl ems, as 
was proposed initially bv fChen~ & Haloer5 (1 9891 1 in the 
form of the ion torus ( an early version of current RIAF 
models; see lR.ees et al.l[l98i3) . 

While spectral models of RIAFs have successfully 
accounted for the spectral energy distrib utions ( here¬ 
after SEDs) of a handful of LINERs (lL a . sota, et aTl 
T996t IQuataert et all 119991 IPtak et all 12004 but see 
Yuan et alJ l2002t iHerrnstein et all |2005 in the case of 
NGC 4258), RIAFs have yet to be studied in a double- 
peaker LLAGN. The reason for this is that there was 
no object until now which had both good observations 
of the nuclear SED and displayed disklike broad emis¬ 
sion lines. We have now such observations for NGC 1097, 
which have allowed us to investigate the structure of the 
accretion flow in the nucleus by comparing the observed 
broadband SED with predictions from theoretical mod¬ 
els in which the plasma at small radii is in the form of 
a RIAF and the outer parts of the flow are a standard 
thin disk. We explore the conditions under which the 
parameters of the thin disk inferred from the Ha profile 
model presented by SB03 are consistent with the mod¬ 
elling of the continuum. In § [21 we describe the mul- 
tiwavelength data comprising the SED, including new 
X-ray observations as well as previously published data 
in the optical/ultraviolet, infrared and radio wavebands. 
In § [3] we derive several properties of the SED, such as 
the degree of radio-loudness and the bolometric and X- 
ray luminosities, and we compare the SED of NGC 1097 
with the typical SEDs of quasars and LLAGNs. In § Q] 
we discuss the accretion flow model and compare its pre¬ 


dictions with the observations. We discuss our results in 
§ Inland present our conclusions in § El Throughout this 
paper we assume a distance 1 to NGC 1097 of 14.5 Mpc, 
thus 1" = 70 pc. In an appendix, we present a spec¬ 
tral analysis of the ~ 1 kpc starforming ring and of an 
ultraluminous compact X-ray source located outside the 
ring. 

2. OBSERVATIONS 

In order to construct the SED of the nucleus of 
NGC 1097, we have used X-ray observations from the 
Chandra X-ray Observatory and near-UV and optical ob¬ 
servations from the Hubble Space Telescope (HST). We 
have also collected infrared and radio fluxes from the lit¬ 
erature. The sections below describe the data in detail. 

2.1. X-RAY DATA 

NGC 1097 was observed with the backside illuminated 
ACIS-S3 CCD aboard t he C handra X-ray Observatory 
on June 28, 2001. Table EH gives a log of the Chandra 
observations. 

Figure ^ shows the image of the central region of the 
galaxy in the 0.5-8.0 keV band, which illustrates the ex¬ 
tended emission morphology. The nucleus (black spot 
in the center) and circumnuclear region (including the 
star-forming ring at ~ 9") are clearly seen in the im¬ 
age. An ultraluminous compact X-ray source (ULX) is 
observed at 27'/6 NNE from the nucleus. The properties 
of the star-forming ring and ULX are discussed in the 
Appendix. 
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Fig. 1.— X-ray image of the nucleus of NGC 1097 (0.5 - 8 keV). 
The angular size of the region shown corresponds to 90 /r x 90 /; and 
the linear size is 6.3 kpc. The ULX coordinates are 02h 46m 19.81s 
-30° 16’ 0278 (J2000). 

An X-ray spectrum of the nucleus was extracted from 
a circular region of radius 276. The effect of pile up is too 

1 We note that IStorchi-Bergmann et af ] <2005) assumed a dis¬ 
tance to NGC 1097 of 17 Mpc. We choose the distance 14.5 Mpc 
in order to com pare directly our results for the SED with those of 
IHo et all J200a) . who adopt the latter distance to NGC 1097. Fur¬ 
thermore, this difference in the assumed distance does not change 
our conclusions and its effect is within flux calibration uncertain¬ 
ties. 
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large in the first observation to measure reliable spectral 
shape and flux. The second observation is also affected 
by pile up, but the effect was relatively small because 
the detector was operated in 1/8 sub-array mode. We 
accounted for pile-up by applying the pile-up model im¬ 
plemented in the XSPEC 2 spectral-fitting package. 

Spectral fits were performed to the X-ray spectrum, 
using a y 2 minimization technique. Acceptable fits were 
obtained using a power-law model (I® a E~ T , where 
Ie is a photon flux in units of photons s _1 cm' 2 keV 1 
and T is the photon index) modified by photoelectric 
absorption along the line of sight. Figure Q shows the 
observed spectrum of the nucleus with the best-fit model 
superposed. The model parameters for this fit are listed 
in Table IA2I where /Vh is the hydrogen column den¬ 
sity and J^-iokeV is the observed flux in the 2-10 keV 
band corrected for absorption. The value of the photon 
index is consistent with the typ ical values observed in 
LLAGNs (T = 1.6 — 2.0; s ee, e.g. lTerashima et alll2()()2t 
lT si.as hima.fc Wilsonll2003I. The X- ray spectrum of the 
nucleus does not have a high enough signal-to-noise ratio 
to constrain spectral features such as Fe K emission. 
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Fig. 2.— X-ray spectrum of NGC 1097. The upper panel shows 
the data and the best-fit model, and the lower panel shows the 
residuals of the fit. 

The ASCA obse rvation ilvomoto et alJ 119961 
iTerashima et alj [2002) includes all of the nucleus, 
the starforming ring and the ULX due to its poor 
angular resolution. The hard component detected 
with ASCA is most likely dominated by the nucleus, 
because the number of 2-8 keV photons from a few faint 
off-nuclear sources and the weak hard emission from 
the starforming ring is about 10% that of the nucleus, 
thus the ASCA spectral fits are reliable measures of 
the spectral shape of the nucleus. The hard band flux 
obtained with ASCA (2.07 x 10 -12 erg s -1 cm -2 in the 
2-10 keV band or 1.71 x 10~ 12 erg s -1 cm -2 in the 
2-8 keV band, corrected for absorption) is consistent 
with the Chandra nuclear flux if the ASCA calibration 
uncertainty (10 %) and a small contribution from the 
off-nuclear hard X-ray emission are taken into account. 

2 http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec 


2.2. OPTICAL/UV DATA 

The nucleus of NGC 1097 was observed with the Space 
Telescope Imaging Spectrograph (STIS) in February 2001 
in the wavelength range 1000 A-l /im, using a 0"2 wide 
and 52" long slit. The details of these obser vations are 
described in iStorchi-Bergmann et alJ (2005) (hereafter 
SB05), who have used these data to infer the presence 
of a starburst within 9 pc from the nucleus. The nuclear 
spectrum used in the present work was extracted using 
a window of 0"2 x 0"6, the latter corresponding to 12 
pixels along the slit. Figure 1 of SB05 shows the nuclear 
spectrum, where the characteristic broad double-peaked 
emission lines can be seen. 

In order to obtain the intrinsic optical-to-UV (OUV) 
nuclear continuum, we have to correct the spectrum for 
the effects of dust extinction. The Galactic extinction 
a ffecting t h e HST spectr um is E(B — V) = 0.027 mag 
llSchlegel. Finkbeiner fc DavisIlOQStl . On the other hand, 
the neutral hydrogen column density measured from the 
X-ray spectral fit is Nn = 2.3 x 10 20 cm -2 (§ 12.111 . Us¬ 
ing_th£j^lati£n_F( ; B — V ) = iV H/(5.8 x 10 21 cm -2 ) mag 
llBohlin. Savage fc Drakd 9978 1 we obtain E(B — V) = 
0.04 mag. Subtracting the contribution of the Galac¬ 
tic reddening above, the resulting internal reddening is 
E(B — V) = 0.013 mag, which is very small. Adopting 
Ry = Ay / E(B — V) =3.1 the total reddening affecting 
the observations is Ay = 0.124 mag, which implies only 
small extinction corrections to the spectrum using the 
Galactic reddening law. 

2.3. INFRARED AND RADIO DATA 

Data in the infrared (IR) and radio wer e collected from 
the literature and are listed in Table f~\Til NGC 1097 has 
a luminou s circunmuclear ri ng of st a r formation of ra - 
dius « 9" (IStorchi-Bergmann. Wilson fc Baldwiri ; 1996 1. 
Most data available in the literature were obtained using 
an aperture which includes this circumnuclear ring, and 
the fluxes so measured can only be considered as upper 
limits to the nuclear emission. Therefore we decided not 
to include such observations in the SED, including those 
made with the IRAS and ISO satellites (apertures in the 
range 24" — l 7 ). The only data points we include in the 
SED are those corresponding to apertures smaller than 
10 ". 

3. THE SPECTRAL ENERGY DISTRIBUTION 

We show in Figure [HI the complete SED obtained by 
combining the Chandra, HST, infrared and radio data 
obtained through apertures smaller than 10". We notice 
that there is a discontinuity between the IR and optical 
data. We believe this is due to additional components 
not included in our HST data (e.g. hot dust), and thus 
these points are treated as upper limits to the nuclear 
emission and are represented by open circles. The filled 
circle is a radio data point observed with a small beam 
(0"25). 

In order to estimate the various quantities which char¬ 
acterize the SED, we must deal with regions of the spec¬ 
trum for which we have no data available. Thus, we 
interpolate the SED with six power-laws of the form 
F v oc v~ a , one for each relevant frequency interval, as 
shown by the long-dashed line in Figure [HI The radio to 
infrared power-law was obtained by connecting the radio 

























4 



Fig. 3. — Spectral energy distribution of the nucleus of NGC 1097 and power-law repre sentations of the SED ( long-dashed line ), as 
described in the text. For comparison, the average SED of LLA GNs (dotted line , |Hq 1I2005H and the average SEDs of radio-loud ( short- 
dashed line) and radio-quiet (dot-dashed line) quasars JeI vi sll994l) are shown. The spectra are normalized at NGC 1097’s nuclear luminosity 
at 2 keV. The open circles denote data possibly contaminated by the circumnuclear starburst ring or host galaxy. 


luminosity value at 8.4 GHz to the average value of the 
normalized LLAGN SED at 2.5 x 10 13 Hz. The infrared 
to optical power-law was obtained by connecting the for¬ 
mer point to the optical luminosity value at 1 /im. The 
OUV power-laws were obtained by fitting the data us¬ 
ing least squares. The power-law between 2.5 x 10 15 Hz 
and 2 keV was obtained by connecting the former point 
to the latter one. The best-fit X-ray power-law is used 
between 2 keV and 100 keV, above 100 keV the SED is 
described by an exponential cutoff. 

3.1. Properties of the SED 

We use the interpolated piece-wise power-law described 
in the previous section to calculate several properties of 
the SED: the radio loudness, bolometric luminosity, X- 
ray luminosity and various useful spectral indices, as fol¬ 
lows. 

Radio loudness - The radio loudness can be quantified 
using two different parameters: an o ptical t o radio ra¬ 
tio d efined as R 0 = R„(6 cm)/L J/ (B) llKellermann et 733 
fTosO i and a radio to X-ray ratio defined as R x = 
^L„(6 cm)/L x , w here L x is the luminos ity in the 2 - 
10 keV band IjTerashima fc Wilsonl 121)1) 3). The bound¬ 
ary between radio-loud and radio-quiet objects corre¬ 
sponds to R 0 = 10 and R x = 3.162 x 10 -5 . NGC 1097 is 
marginally radio-loud according to both criteria ( R 0 = 
11, -R x = 7.6 x 10~ 5 '), much like other LLAGNs pre¬ 
vi ously studied (|Ho| [l99£; IHo fc Pendl20ml lHcll2002al 
ITerashima fe Wilsonl2005 ). IHo et alJ 1120001 have quoted 
a value of R 0 = 19, indicating a somewhat larger degree 
of radio-loudness, but they do not specify what data they 
used to derive this value. 

Bolometric luminosity - The bolometric luminosity is 
AboI = 8.6 x 10 41 erg s _1 . This value is well within 
the range observed for o ther L LAGNs (2.1 x 10 41 — 8.0 x 
10 42 erg s _1 , from Ado pting a black hole 

mass of 1.2 X 10 8 /l/7 Thewis X Eracleondl2005 ). we ob¬ 
tain an Eddington ratio of RBoi/iEdd = 5.7x 10 -5 , where 
^Edd = 1-25 x 10 38 to = 1.5 x 10 46 erg s _1 , (where to 


is t h e bl ack hole mass in units of Mg). For comparison, 
fihl quotes a value for this ratio in the interval 

~ 10 -6 — 10 -3 for LLAGNs. RboI and Tsoi/REdd are 
app roximate l y two times smaller than the values quoted 
bv lHo et~ ( 1200 01. 

X-ray luminosity - The X-ray luminosity in the 2-10 
keV band is L x = 4.4 x 10 4 ° erg s _1 , which is in the 
i nterval of X-ray luminosi ties of LLAGNs studied by 
ITerashima. fc Wilsonl (120081) . 5 x 10 38 — 8 x 10 41 erg s _1 . 
The ratio of L x to narrow Ha luminosity is 300, where the 
narrow Ha luminosity is 1.46 x 10 38 erg s _1 , measured 
from the HST spectrum. This rat io somewhat exceeds 
the range of values estimated by ITerashima & Wilsonl 
<20031) . 0.27 104. The ratio of X-ray luminosity in the 

0.5 - 10 keV band to Lro) i s 0.08 , which is within the 
range of values found bv llloi 119991) for LLAGNs, 0.06 
0.33. 

Spectral indices - The value of the energy index of the 
continuum in the hard X-ray band (2-10 keV) can be 
obtained from the photon index (§ l2~fll using the relation 
a x = T — 1 = 0.64. We measured the value of the power- 
law index a ou of the OUV continuum in § 0 as being 
1.87. The two-point spectral index a ox from 2500 A- 2 
keV is 1.14, which is similar to the mean value found by 

153 (fi999l) for LLAGNs, (a ox ) ~ 0.9. 

3.2. Comparison with SEDs of other AGNs 

Figure^shows the SED of NGC 1097 together with the 
average SED of LLAGNs from iHoi ( 2005 ) and a verage 
SEDs of radio-loud and radio-quiet quasars from lElvisl 
fl99l . All curves are normalized at the 2 keV luminosity 
of NGC 1097. 

In Figure 0 it can be seen that in spite of NGC 1097’s 
SED being more radio-quiet than the average SED of 
LLAGNs, the general shape of NGC 1097’s SED is simi¬ 
lar to this mean SED. This result has been verified quan¬ 
titatively in the previous section through the similarity 
between the values of the spectral indices from the two 
SEDs. It can also be seen in Figure|3]that the continuum 
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of NGC 1097 lacks the canonical UV bump commonly 
present in quasars, which is associated with emission 
fro m a geometrical ly t hin, optically thick accretion disk 
llFrank. King fc Bainelf2002l i. This was already not ed in 
the particular case of NGC 109 7 by Ho et al ■ ™ and 
more generally in LLAGNs bv iHol ( 1999 1. and indicates 
an absence or truncation of the inner thin accretion disk 
in these objects. The latter interpretation will be favored 
in the modelling of the next section. 

4. MODELS FOR THE SED 

The observational characteristics of NGC 1097 de¬ 
scribed in Se ction § 0 (and in general of LLAGNs; 
see iHol xM for a review) can be understood within 
the theoretical framework of RIAFs (see reviews by 
iQua.ta.ert] 1200 ll IN a.ravaiil 120051 early versions of RIAF 
models were called ADAFs or “ion tori”; see NMQ98, 
iKato. Fukue & Mineshigelll998i for reviews'). RIAFs are 
optically thin, geometrically thick accretion flows that 
have radiative efficiencies much less than the canonical 
10% of thin disks and thus naturally generate low lu¬ 
minosities. RIAFs are thought to arise when the mass 
accretion rate is below a critical value Merit ~ O.OlMEdd 
(NMQ98). A low Eddington ratio is a requirement for 
the existence of RIAFs (ABoi/AEdd < 0.01 — 0.1; see, e.g., 
NMQ98) and the low ionizing luminosities produced by 
RIAFs coupled with the typical densities of the narrow 
line region lead naturally to low values of the ionization 
parameter U , which consequently give rise to LINER- 
like spectra llHaloern fc Steine~Hll98.fl: iFerland fc Netzeil 
1983). 

In view of these indications, the model we adopt 
consists of an inner RIAF plus an outer standard 
thin disk, which has been successful in reproduc¬ 
ing the SEDs of different kinds of sources, including 
black hole X-ray binaries, Sgr A*, LLAGN s in ellip¬ 
tical galaxies and LI NERs (for reviews, see IQuataert! 
120011 iNaravanl 120051) . We incorporate in the mod¬ 
elling of the RIAF recent theoretical ad vances. As al¬ 
ready p redict ed in the pioneering papers lNa,ra,van fc Yil 
( 1994 1*9951). and confirmed by recent numerical 
simulations fcton e et alJ llj)99yHawleY_fc^albiijj2 002l: 

llgumenshchev et, al.1120031 iProga &; Begelmanll2Q0a) and 

analytical calculations ll Blandford fc Begelmanl 1 19991 

iNaravan et alJl200fi: iQuataert fc Gruzinovll2000alh only 
a small fraction of the gas supplied at large distances 
from the RIAF actually falls on to the black hole, the 
rest of the gas is either ejected in an outflow (advection- 
dominated inflow/outflow solution [“ADIOS”]) or is pre¬ 
vented from accreting by convective motions (convection- 
dominated accretion flow [CDAF]). 

Given the theoretical uncertainties on the struc- 

tu re and micr o physi cs_of the flow, we follow 

iBlandford Kr. Regelman! (1 9991) and assume that the mass 
accretion rate varies with radius as M(r) = Mq (r/rt T ) p , 
where Mo is the accretion rate at the transition radius 
rtr between the inner RIAF and the outer thin disk. 
The exponent p and the fraction of the turbulent energy 
that heats the electrons, 5, are free parameters in our 
model. This assumptio n has been recently confir med by 
MHD simulations (e.g., IProga fc Begelma nll2003D. 

We assume m = 1.2 x 1O 8 M 0 llLewis fc Eracleousl 
l2005lh r out = 10 5 , where r ou t is the outer radius of the 


thin disk, and i = 34° (SB03) is the inclination of the 
axis of the disk to the line of sight. We set the plasma 
parameter (ratio of gas to total pressure) to the value 
0 = 0.9 (for motivation, see lQuataert fo Naravanll999bfl 
and we determine the value of the viscosity parameter a 
from P usin g the equation a ss (3/2')(l — 0)/(3 —20) « 0.1 
(NMQ98, iHawlev. Gammie & Balbusl 11996 ). We ad¬ 
just the accretion rate to reproduce the observed X- 
ray flux, and hereafter we use the accretion rate in Ed¬ 
dington units 771 (rh = M/M^dd)- We solve the radia¬ 
tive transfer and hydrodynamical equations of the RIAF 
self- consistent ly, obtaining the e mergent spectrum (see 
[Yuan. Quataert fo Naravanl2003l for details). The radia¬ 
tive processes considered in the RIAF are synchrotron 
and bremsstrahlung emission, and Comptonization of 
these photons and the soft photons from the outer thin 
disk. The thin disk is assumed to emit locally as a black- 
body and reprocesses X-ray photons from the RIAF. 

The value of the transition radius was obtained from 
fitting the double-peaked Ha line profile with disk models 
as w 225 Rs (SB03). On the other hand, there is also an 
independent constraint on r tr from fitting the SED with 
the accretion flow model, i.e. the value of r tr must be 
suitable to fit the OUV spectrum. We set the value of 
r tr to be that obtained from the line profile fitting, i.e. 
r tr = 225. 

In Figure 0] we show the predicted spectrum of the 
accretion flow for rho = 6.4 x 10 -3 (Mo in Eddington 
units), p = 0.8 and S = 0.1; the solid line shows the RIAF 
contribution and the dashed line shows the emission of 
the truncated thin disk. It is encouraging that we obtain 
a good fit of the optical and X-ray portion of the SED 
using the value of r tr = 225 derived from the line profile 
model, i.e. our models for the SED and the line profile 
are consistent with each other! _ 

As shown bv IQuataert fc Naravanl (I 1999afb when com¬ 
paring spectral models of RIAFs with observations, there 
are degeneracies between the mass-loss rate in the wind 
(parameter p) and parameters describing the micro¬ 
physics of the RIAF, particularly <5, so the fits to the data 
presented here are not unique; instead they should be 
treated as representative models. For example, we can fit 
the SED equally well with an ADAF model, where p = 0 
and 6 = 0.01. It is worth noting t hat there are theoreti¬ 
cal uncertainties on the value of S (IQuataert fc Gruzinovl 
EH) , and that as the accretion rate increases, it is ex¬ 
pected that outflows/convection become less well devel¬ 
oped, because the Bernoulli parameter decreases with 
increasing rh. 

The accretion flow model severely underpre¬ 
dicts the radio emission and underpredicts to a 
small degree the higher frequency UV data. Un¬ 
derprediction of radio data is common in RIAF 
models (e.g., Quataer^et^L I IT999t IPtTk et all l20?il 

jlJlvestad fc Hoi 120011: lAnderson etahl 12004 IWu fe Cad 

l2005|) , and the radio emission would be enhanced if 
the s ource has a jet ( e.g.. lYuan. .^arkofffe Rdckd 
2002 : [Yuan~et~ . Boo)! lYuan. Cui fc Naravanl 120051 
hereafter YCN05) or if the RIAF has nonthermal 
electr ons dMah adevanl Il999t lOzel. Psa itis V Naravanl 
I2000t lYuan. Quataert fc Naravanl 120031) . To test the 
possibility that the observed radio emission comes 
from a jet, we use the model for the jet emission of 
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Fig. 4.— Models of the RIAF ( short-dashed, line), thin disk ( long-dashed line), jet (dotted line) and obscured starburst ( dot-dashed line) 
compared to the nuclear SED of NGC 1097. The sum of all components is also shown ( solid line). The thin disk is truncated at rtr = 225, 
inside of which there is a RIAF; the accretion rate decreases inwards according to m(r) = 6.4 X 10“ 3 (r/rtr) 0 ' 8 (see text). The starburst 
model includes the Fe II emission-line. 


YCN05, which adopts the internal shock scenario widely 
employed in the interpretation of gamma-ray burst 
afterglows. Since the only constraint on the jet emission 
is the observed radio flux at 8.4 GHz, we assume that 
the most basic parameters have the same values as in 
YCN05. Namely the half-opening angle is <j> = 0.1, the 
energy density of accelerated electrons is e e = 0.06, and 
the amplified magnetic held energy density is e_e = 0.02. 
The bulk Lorentz factor of the jet is Tj = 10, the 
viewing angle is 34° (same as the inclination angle of the 
thin accretion disk inferred by SB03) and the mass-loss 
rate in the jet is M- ]et = 7 x 10~ 7 MEdd, which is about 
0.5% of the mass accretion rate of the ADIOS at 5 Rs- 
The resulting jet spectrum is shown in Figure 0] as the 
dotted line. We see that this jet model accounts for the 
available high-resolution radio data. The radio band 
is dominated by the jet, and the jet emission becomes 
less important in the rest of the SED compared to the 
accretion flow emission. We note that as we have so few 
constraints on the jet emission, the jet parameters are 
not as well constrained as those of the accretion flow 
models, although the results are not very dependent on 
the values of the jet parameters. 

The discrepancy between the model and the observed 
UV flux can be resolved by invoking the contribution of 
a compact starburst to the observed UV spectrum, as 
found by SB05. The starburst has a mass of approxi¬ 
mately 10 6 solar masses, an age of at most a few x 10 6 
yrs and is obscured by Ay ~ 3 mag. The signatures of 
the starburst include a number of absorption features in 
the 1000-1600 A spectral region and the continuum con¬ 
tribution illustrated in Figure 0 as the dot-dashed line. 
Included in the theoretical SED is an Fe II template 
broadened according to the width of the Balmer lines, 
which accounts for the observed bump around 2500 A 
(see SB05 for details). This UV Fe II emission presum¬ 
ably comes from the thin disk. 

We note that while the starburst is reddened by Ay = 


3 mag, the X-ray data, the double-peaked lines (including 
the UV Mg II lines) and our successful modelling suggest 
that the AGN is unobscured, supporting the suggestion 
by SB05 that the starburst could be in a circumnuclear 
dusty structure, which does not intercept our line of sight 
to the nucleus. 

We also notice (see FigureQ} that our best model some¬ 
what overpredicts the luminosity of the UV end of the 
HST spectrum. This part of the spectrum is very sen¬ 
sitive to the amount of dust extinction. If the value of 
Ay is slightly larger, then the intrinsic UV flux would 
be higher than the adopted one and would agree bet- 
ter with the model. Interestingly enough, iPrieto et alJ 
l>2005l) measured Ay ~ 1 within the central ~ 60 — 40 
pc region, a higher extinction than the value adopted in 
this work. We can obtain a small improvement in the fit 
of the UV end varying the parameters rho and p of the 
RIAF model, but in that case the quality of the fit at the 
red end of the HST spectrum will worsen. 

5. DISCUSSION 

The adopted transition radius (rt r ~ 225) is similar to 
the typical values found in available RIAF + thin disk 
mod els of LINER SEDs (NGC 4579 and M81: r tI « 
i nn. IQua.ta.ert et al.lll999t NGC 4258 : r tr « 10 - inn, 
iGammie . Narava n fc Blandfordlll99?l NGC 3998: r tr — 
3nn. tPtaJc.et alll2004D . 

We tried to fit the OUV data with a model of a thin 
disk extending down to the ISCO at 3 Rs- Although 
a model with m = 3x in -5 fits the UV data, it cannot 
account for lower frequency optical data and it would im¬ 
ply that ss 59% of Lboi is due to X-rays (possibly emitted 
by a corona ab ove the disk) , something very unusual for 
an AGN (e.g.. iHaardt, V Ma,ra,schi|lll993h . This fit shows 
that truncation of the thin disk is definitely needed in 
order to obtain a consistent picture for the nature of the 
low optica l flux, as found in similar models of oth er LIN- 
E Rs (e.g- .lQuataert et alJ ll999t iPtak et ai~ll2004D . 

I Yuan fc Naravanl ~ il2004l) (hereafter YN04) assembled 
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a plot showing the dependence of I/boi/^Edd on rt r for 
many accreting black hole sources fitted with RIAF mod¬ 
els, including for example Sgr A*, X-ray binaries and 
LLAGNs (Figure 3 of YN04). This plot relates two fun¬ 
damental properties of accreting black holes and provides 
a way of putting in context these sources. Interestingly 
enough, when the corresponding values of NGC 1097 are 
inserted in this plot, we find that NGC 1097 falls within 
the relation found by YN04, and is located very close to 
other LINERs previously m odeled with RIAFs (M81 and 
NGC 4579. lQuataert et alJlll99 ffl. 

Using the estimated accretion rate from the accretion 
flow model and the black hole mass, we can make an es¬ 
timate of the gravitational power output of the thin disk 
as Wdisk = 2.1 x 10 42 erg s~ 4 , which is higher than the 
Ha luminosity by two orders of magnitude. Although 
this calculation suggests that in the particular case of 
NGC 1097 the energy budget of the thin disk is such that 
the disk could power locally its own emission lines, the 
direct conversion of gravitational energy into line emis¬ 
sion is improbable, because most of the binding energy of 
the thin disk is radiated in the form of a weakly ionizing 
IR/optical bump which is not capable of photoionizing 
the atmosphere of the thin disk and therefore cannot 
contribute to the formation of the broad double-peaked 
emission lines. Furthermore, there are theoretical un¬ 
certainties in the accretion flow model which affect the 
determination of lin. 

Current observations do not allow us to distinguish 
between the different RIAF models, although there are 
stro ng th eoretic al reasons favoring ADIOS / CDAF (e.g., 
iQuataertHhOO lI over to ADAF models. There are some 
future observations that may help us assess the im¬ 
portance of outflows in RIAFs. For instance, the ab¬ 
sence of linear polarization in the high-frequency ra¬ 
dio spectrum of NGC 1097 would argue in favor of 
ADIOS/C DAF models, even in the presence of a jet 
llQuataert fc Gruzinovi2000bll . This would be attainable 
with high-frequency VLBI interferometry. More radio 
observations to constrain the shape of the radio spec¬ 
trum are also needed to verify whether it is consistent 
with the predictions of the assumed jet model. 

6. SUMMARY AND CONCLUSIONS 

We have reproduced the nuclear optical to X-ray SED 
of NGC 1097, composed of HST and Chandra observa¬ 
tions, using a model in which the accretion flow consists 
of an optically thin, geometrically thick RIAF for radii 
smaller than the transition radius rt r = 225 Rg and an 
optically thick, geometrically thin disk for larger radii. 
The value of rt r is the same as the inner disk radius 
obtained from fitting the double-peaked Ha line profile 
(SB03) and the value of the black hole mass is that in¬ 
ferred from the stellar velocity dispersion of the host 
galaxy. Thus, an important result of this work is the 
consistency between the result of fitting the SED with 
the accretion flow model and the double-peaked Balmer 
line with the thin disk model, i.e. we obtain acceptable 
fits of these two different data sets with the same transi¬ 
tion radius. 

The accretion flow model underpredicts slightly the 
UV continuum, but the added contribution of a young 
and obscured starburst of mass ~ 10 6 M 0 (SB05) can 
account for the excess UV emission. The accretion model 


also underpredicts the radio emission, suggesting the 
need for a jet component. We have been able to re¬ 
produce the radio data with a jet model employing the 
internal shock scenario (YCN05), although the jet pa¬ 
rameters are not as well constrained as those of the ac¬ 
cretion flow models, because we have few constraints on 
the jet emission. 

The results of this work support the proposal that the 
basic components of the nuclei of LLA GNs are a RIAF, 
a truncated thin disk and a jet (e.g., fH5t2005ll . with 
the possibly important role of starbursts presumably ob¬ 
scured by the circumnuclear torus predicted in the uni¬ 
fied model of AGNs. 
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APPENDIX 


X-RAY SPECTRA OF CIRCUMNUCLEAR STAR-FORMING REGION AND ULTRALUMINOUS 

COMPACT X-RAY SOURCE 

X-ray spectra of the circumnuclear star-forming region at ~ 9 ,; from the nucleus 3 and the ULX were taken with the 
Chandra X-ray Observatory. These were analyzed as described in the main text, with the following results. 

Figure fXTI shows the spectrum of the circumnuclear star-forming region. Since the emission is diffuse, the effect of pile 
up is negligible for both observations. We performed spectral fits to the combined spectrum of the two observations. 
We applied a MEKAL plasma model to the spectrum, but the resulting fit was unacceptable, with significant positive 
residuals remaining above 2 keV. For this fit, we obtained Ah = 1.2 x 10 21 cnr 2 , kT = 0.63 keV, and an abundance 
of 0.099 times solar (with y 2 = 114.0 for 62 dof). 

Since the simple MEKAL model fit was poor, we added a hard component to the model, represented by a thermal 
bremsstrahlung model with kT = 7 k eV. This mode l has b een successfully applied in the case of superposition of 
emission from low-mass X-ray binaries llMakishima et al!l989f) and r esult s in a better fit to the spectrum, as illustrated 
in Figure EH The parameters for the spectral fit are listed in Table E3 

The ULX h as a s pectral shape and flux in the two observations which are consistent with each other within statistical 
errors. Figure ^21 shows the spectrum of the ULX. Spectral fits were performed to the combined spectrum of the two 
observations using a maximum-likelihood method. A power law model gives a better fit than a m ulti-color disk 
blackbody model (MCD). The resulting parameters for the power law model are listed in Table E21 assuming the 
object is associated with NGC 1097. 



0.5 1 2 5 10 

Energy (keV) 


Fig. Al.— X-ray spectrum of the circumnuclear starforming region of NGC 1097. The lower dot-dashed curve represents the thermal 
bremsstrahlung component, the upper dashed one represents the MEKAL plasma component and the solid line is the composite model. 
The lower panel shows the residuals of the fit of the composite model. 
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Fig. A2.— X-ray spectrum of the ultraluminous X-ray source (ULX) of NGC 1097. The solid line represents the power-law model. 


3 Note that this is not the nuclear starburst found by SB05 
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TABLE A1 

X-RAY OBSERVATION LOG. 


Observation ID 

Exposure (s) a 

CCD Mode 

Nuclear Counts' 3 

Count Rate 
per Second 

1611 

5343 

Full array 

755 

0.14 

2339 

5442 

1/8 sub-array 

1690 

0.31 


a Effective exposure time after data screening. 

b Counts obtained in the nuclear region (r = 2 / . / 6), not corrected for pile up. The number 
of counts in the first observation is small because of significant pile up. 


TABLE A2 

Power-law fits to X-ray spectra 



_/V H 


-F2-IO keV 


Component 

(xl0 2U cm -2 ) 

r 

(10 _i2 erg s -1 cm -2 ) 

X 2 /dof 

Nucleus a 

o q+2.8 
Z.O—1 7 

1 64 +0 ' 13 
0.07 

1.73 

60.5/70 

Star-forming ring b 

7 7+3. i 

'•* -3.0 

- 

0.151 c 

63.4/61 

ULX 

q+3.8 

U -0 

1 1515+0.37 
1,oo -0.18 

0.14 

11.9/13 


a Corrected for pileup. The normalization constant is 3.88 ^q '^ 9 X 

10 -4 photons keV —1 cm -2 s —1 at 1 keV (after correction for absorption). The errors 
quoted represent the 90% confidence level for one parameter of interest (A % 2 = 2.7). 

b MEKAL plasma + thermal bremsstrahlung model. The best-fit parameters for the 
MEKAL component are kT = 0.61^qq 4 keV and an abundance of O.I 8 +QQ 5 solar. kT = 
7 keV is assumed for the thermal bremsstrahlung component. 

c Bremsstrahlung component only. The MEKAL component contributes an observed 
flux in the 0.5 — 4 keV band of 5.04 X 10 “ 13 erg s - 1 cm -2 (corrected for absorption). 


TABLE A3 

Radio and infrared data for the nucleus of NGC 1097. 


v (Hz) 

S v (Jy) 

Uncertainty (Jy) 

vL„ (erg s 4 ) 

Resolution (”) 

Reference 

8.4 x 10 9 

0.0031 

4.51 x 10~ 4 

1.24 X 10 36 

0.25 

1 

1.5 x 10 13 

0.240 

0.073 

8.89 x 10 41 

5 

2 

2.94 x 10 13 

0.065 

0.009 

4.67 x 10 41 

5 

2 

8.47 x 10 13 

0.0731 

- 

8.71 x 10 41 

6 

3 

1.35 x 10 14 

0.0409 

- 

2.48 x 10 42 

6 

3 

1.39 x 10 14 

0.0802 

- 

2.80 x 10 42 

0.18 x 0.15 a 

4 

1.87 x 10 14 

0.0878 

- 

4.14 x 10 42 

6 

3 

1.80 x 10 14 

0.049 

- 

2.22 x 10 42 

0.20 x 0.18 a 

4 

2.40 x 10 14 

0.0604 

- 

3.64 x 10 42 

6 

3 

2.43 x 10 14 

0.023 

- 

1.41 x 10 42 

0.20 x 0.19 a 

4 


Refer ences. — (1) |ThearL et alj|200Cl (2) ITelesco &; G atlcv 19^ (3) IGlass &; Moorwood 
119851 Ml IPrieto et all20r)Hl 

a The quoted aperture is the FWHM of the achieved spatial resolution in the respective band. 





















